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Abstract 

The effects of radiation reaction force on laser driven auto-resonant particle acceleration scheme 
are studied using Landau-Lifshitz equation of motion. These studies are carried out for both linear 
as well as circularly polarized laser fields in the presence of static axial magnetic field. From the 
parametric study, a radiation reaction dominated region has been identified in which the particle 
dynamics is greatly effected by this force. In the radiation reaction dominated region the two 
significant effects on particle dynamics are seen viz., (1) saturation in energy gain by the initially 
resonant particle, (2) net energy gain by a initially non-resonant particle which is caused due to 
resonance broadening. It has been further shown that with the optimum choice of parameters 
this scheme can be efficiently used to produce electrons with energies in the range of hundreds 
of TeV. The quantum corrections to the Landu-Lifshitz equation of motion have also been taken 
into account. The difference in the energy gain estimates of the particle by the quantum corrected 
and classical Landu-Lifshitz equation are found to be insignificant for the present day as well as 
upcoming laser facilities. 
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I. INTRODUCTION 


The rejuvenated interest in the Auto-Resonant scheme which in principle can cause large 
unbound acceleration of particle is due to the concurrent developments in the ultra high 
magnetic field generation experiments and laser technology. The advancement in these ar¬ 
eas have lead to the availability of non-destructive magnetic helds of the order of (~ 100 
Tesla) P-[3] together with lasers having intensities of the order of lO^^fU/cm^ |1] which are 
expected to increase even further ISHH. The scheme of Auto-Resonant particle acceleration 
was discovered by Kolomenskil and Lebedev [HI [9] and, independently, by Davydovskil [10] . 
The previous studies [ [TTfITF] ] on the topic have not considered the effect of radiation reaction 
in this scheme. The radiation reaction (or radiative friction) is a recoil force acting on the 
accelerating (retarding) particle as a result of emission of radiation and can in turn influence 
the motion of the particle. In the absence of radiation reaction it is known that there is no 
net transfer of energy to the particle interacting with the laser and non-resonant magnetic 
held. However it has been recently shown from exact analytical results [28] that the inclusion 
of radiation reaction which has an amplitude dependent coefficient of friction results in the 
net transfer of energy to the particle interacting with laser held. In this context, the theo¬ 
retical understanding of the ehect of radiation reaction force on the Auto-resonant scheme 
can be of central importance in the optimum use of these laboratory generated ultra high 
magnetic helds and lasers for the experimental realization of the scheme. 

In the Auto-Resonant acceleration scheme, the particle is subjected to the combined 
held of laser and static axial magnetic held. The particles are accelerated as a result of self- 
sustained resonance between the particle and the laser held due to which the particle remains 
phase locked with the laser held. The resonant phase locking condition can be expressed as, 
T — — = 0, where f2c(= ^^) is the cyclotron frequency normalized to laser frequency cn, 

^Px/md is the normalized longitudinal component of particle four momentum and T' is the 
relativistic factor. The mechanism for particle acceleration by auto-resonance scheme can 
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be briefly described in a following way: the particle initially at rest and satisfying the initial 
resonance condition {u = ffc) is accelerated along the electric held of laser. The increase 
in particle energy along the electric held component in the transverse direction results in 
relativistic mass increase which in turn lowers the cyclotron frequency of the particle. The 
transverse velocity gain results in the action by the magnetic held component of the wave, 
which pushes the particle along the direction of the propagation of the wave. The relativistic 
velocity acquired by the particle along the longitudinal direction results in a Doppler shift to 
a lower frequency of the wave as “seen” by the particle. In this case, the Doppler shifted wave 
frequency to the lower frequency equals to the relativistically modihed cyclotron frequency, 
and the particle remains “synchronously” in cyclotron-resonance condition. This result in 
continuous increase in particle energy and momentum along the direction of propagation of 
the wave. 

The present work takes into account the effect of radiation reaction force (or radiative 
friction) on the dynamics of a particle in the combined held of intense laser and static axial 
magnetic held using the Landau-Lifshitz (LL) equation |37j. The choice of LL equation over 
the Lorentz-Abraham-Dirac (LAD) equation in the study of particle dynamics is because of 
the reason that it circumvents the problems such as violation of principle of causality and 
existence of nonphysical runway solutions associated with the LAD equation[3H]. The LL 
equation can be perturbatively derived from LAD equation [37j and thus is equivalent to it 
in hrst order. This equation has also been used in the recent laser matter interaction studies 
[28H36] to highlight the ehect of radiation reaction force on particle dynamics for the present 
day as well for upcoming laser facilities. 

In the present work radiative friction dominated parametric regime has been explored to 
study the dynamics of initially resonant as well as non-resonant particle interacting with the 
combined held of laser and static axial magnetic held. In this regime, at hrst the energy gain 
of the initially resonant particle is estimated at which the radiation friction ehects on the 
particle dynamics become prominent for the present day as well as upcoming laser facilities. 


3 


For the initially non-resonant particle it is well known that in the absence of radiative friction 
there is no net transfer of energy at the end of each gyration. In the present study the effect 
of de-phasing caused by amplitude dependent radiative friction on the initially non-resonant 
particle from the laser held has also been investigated in the presence static axial magnetic 
held. A detailed comparative study has been carried out on the peak energy gain of the 
particle by linearly and circularly polarized lasers. Further the ehect of quantum correction 
[ini m] in the LL equation of motion has also been considered for estimating the energy gain 
by the particle. The optimum conditions have been predicted for the large peak energy gain 
by the particle. 

The organization of the paper is as follows: In section (II), mathematical frame work for 
studying the particle dynamics in the combined held of laser and static axial magnetic held 
is presented. In section (III), the Landau Lifshitz equation of motion is numerically solved 
using R.K(Runge-Kutta) method for Linear and Circular polarization of the laser. Section 
(IV), contains the summary and discussion of the results. 

II. LANDAU-LIFSHITZ EQUATION FOR PARTICLE MOTION IN THE 
COMBINED FIELD OF LASER AND STATIC AXIAL MAGNETIC FIELD 

The ehect of radiation damping on relativistic particle dynamics is described by the 
Lorentz-Abraham-Dirac (LAD) equation. 


du* e ,■ 

me-— = -F Uk + g 


( 1 ) 


dr c 



( 2 ) 


Here u* = (F,P/mc) is the four velocity, = diA^ — d^Ai is the EM held tensor with 


Ai = ($, A) being the four potential. The radiation friction force in the LAD equation 
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[371 [38] is given by, 


9^ = 


2e^ ,(fu^ 


3c dr^ 




d^Uk 

dr"^ 


(3) 


The LL radiation friction force is derived as a perturbation to the equations of motion, this 
approximation is valid provided there exists a frame of reference, where the LL radiation 
friction force is small compared to the Lorentz force. Using the equation of motion 1371 
(neglecting in the hrst approximation), ^ and are evaluated as 


dn* 

dr 

d^u'' 

dr'^ 


nc 


-.F^^Uk 


p f) p‘^ 

-^^UkU^ + -^r’^Fkju^ 


(4) 

(5) 


Substituting the above expression in the Eq [3] the radiation friction force in the LL form 
stands as 


9 = 


2e^ 

3mc^ dx^ 
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Using the above expression for the LL equation [37| of motion is given by 


dti® 


mc~— = -F'-^Uk + 
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-UkU 
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{FkiF){F^^u^)u^ ( 6 ) 


The 3-D spatial part of the radiative friction corresponding to additional terms in the 
momentum equation in the LL equation is 


/ = 


2e^ 


3mc^ 


(1 - + itV)B + iff X (| + , 7 .V)ff} 


\^E X H +-^H X {H X v) + ^U(U.U)| 




2e^ 


(7) 


X - , +-V X Hf - \{E.vf 

3m^c^{l-%) '' " ^ ^ 

In the following section the numerical calculations are carried out using all the three terms in 
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the radiation reaction force given by above equation. However from the order comparison as 
well as in the numerical simulations the third term has been found to be the most dominant 
term. 

In the absence of radiative friction the inclusion of static axial magnetic held in Eq. [1] 
couples the particle dynamics along the transverse spatial components and decoupling of 
the motion yields driven harmonic oscillator equations along these components [251 Wf\ . The 
oscillators are driven by externally applied held of laser together with the static axial mag¬ 
netic helds. The inclusion of radiative friction causes the dampening of these oscillators and 
non-linear coefficient of friction depends upon the amplitude of the applied held. 

In the present work, the quantum electrodynamics weakening of the radiation friction 
HQHH] has also been taken into account by multiplying the radiative friction part of LL equa¬ 
tion with a form-factor.The form-Factor is given by Ge{Xe)^ where Xe{= ^ ~) 

ativistic and gauge invariant parameter characterizing the probability of the gamma-photon 

3 3 

emission by the electron with four momentum and Es{= is the ‘Sauter-Schwinger' 

held. The 3D spatial part the parameter Xe is given by 

X. = (8) 

For numerical calculations [H] the corresponding form-factor is given by 

“ (1 + 8.93X. + 2.41x2)2/3 

The above equations can be simplihed by introducing a constant tq, with dimension of time 
and is given by Tq = For electron/positron it has a numerical value Tq = 6.24 x 

The equation of motion is numerically integrated by expressing the variables in dimen¬ 
sionless form by using the following normalizations: R{x, y, z) —)■ kR,t —)■ ut, P —)■ T —)■ 

—>■ —)■ -^,A —>■ -24^, He —t -2^. The vector potential of an elliptically po- 

' m-Lur' m.ruj' ^ ni.ru) jt j. 
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larized laser traveling along x direction in the presence of a constant homogeneous axial 
magnetic held is given by, 

A = |^v4o(l - y + 5 (10) 

in the above expression k dehnes the polarization state of the laser for linear polarization 
(k = 0, ±1) and the circular polarization of the laser is dehned by (k = The other 

terms as specihed above are: the phase of the laser is given by ^ = {ut — kx), a(^) = (Sin(^)) 
and o'i(.^)(= Cos(^)) are the oscillatory parts, Aq is the peak laser amplitude and Bq is the 
magnitude of the external magnetic held. The equation of motion has been integrated in 
terms of co-ordinate 

III. RESULTS OF NUMERICAL SOLUTION OF THE LL EQUATION OF MO¬ 
TION 

The dynamics of charged particle in the combined held of laser and static axial magnetic 
held is studied taking radiation reaction force into account. The dynamics of the particle is 
obtained by numerically solving the LL equation of motion using Runge-Kutta integration 
scheme. The particle dynamics is studied for both linearly and circularly polarized laser. The 
numerical results have been validated by comparing them with the exact analytical results 
[25] for particle dynamics in the linear and circularly polarisied laser in the absence of axial 
magnetic held. The results obtained by solving the LL equation are compared with that of 
the known results of Lorentz force equation. In the following study for the sake of convenience 
the previously expressed resonance condition has been re-arranged and expressed in terms 
of parameter ‘r(= ^)^ which is the ratio of cyclotron frequency to the product of laser 
frequency ‘u' with a variable ‘A(= T — PxY- The resonant interaction dehned by ‘r' = 1. 

The conhguration and momentum space trajectories of the initially resonant particle in 
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the absence as well as presence of the radiative friction in the held of circularly polarized 
laser are shown in Fig-1. From the hgure it can be seen, that the radiative friction is not 
prominent in the initial part of interaction. The resonant interaction accelerates the particle 
along the direction of propagation of laser which results in the increase in its orbit size for 
both momentum and conhguration space trajectories. On comparing the momentum space 
particle trajectory in sub-plot (A) and (B), it is found that the effect of radiation reaction 
force on the particle dynamics becomes pronounced with the increase in particle momentum. 
The prominence of radiative friction causes a marked difference between the momentum as 
well as conhguration space trajectories which are shown in sub-plots (A) and (C) respectively. 
From the hgure it can be inferred that the presence of radiative friction in turn pushes the 
particle out of resonance which results in the decrease of the particle momentum. The 
decrease in particle momentum results in the shrinking of its corresponding conhguration 
space trajectory in the transverse plane. 

The comparison of the ehect of inclusion of radiative friction on the initially non-resonant 
particle trajectories (i.e r = .7) in the momentum and conhguration space to that of its 
absence is shown in Fig-2. From sub-plots A and B, it is evident that in contrast to the 
case when the radiation friction in neglected and there is no net transfer of momentum to 
the particle, the inclusion of the radiative friction results in the gradual de-phasing of the 
particle from the applied laser and magnetic held. The de-phasing of the particle results in 
net gradual increase in particle momentum which in turn pushes the particle to the cyclotron 
resonance and resulting in large increase in longitudinal particle momentum. The change 
(increase and decrease) in particle momentum results in the corresponding variation in the 
size of particle orbit in both momentum and conhguration space trajectories. The resonant 
particle as in the earlier case is then gradually pushed out of the resonance. 

From the both trajectories it can be observed that the radiative friction causes the de¬ 
phasing of the particle from the laser held resulting in gradual increase in the particle mo¬ 
mentum and increase in its orbit size. The increase in longitudinal momentum causes a 



gradual Doppler down shift in frequency of the laser as ‘seen' by the particle and thus driv¬ 
ing it in cyclotron resonance with the applied helds. The resonant interaction results in rapid 
increase in particle momentum. The continuous evolution of resonance parameter beyond 
resonant condition given by ‘r' = 1 drives the particle out of resonance with applied helds 
resulting in fall of particle momentum and its eventual saturation. 

In the Fig-3, the plot of the energy gain as a function of null co-ordinate ‘^(= t — x)' is 
shown for the initially resonant as well non-resonant particle with the applied laser and static 
axial magnetic held. From the sub plots A and B of the hgure, it can be inferred that the 
peak energy gain in the case of initially resonant particle for a given laser intensity is larger 
in the case of circularly polarized laser than in comparison to the linear polarization. The 
radiative friction assisted de-phasing of particle from the laser held results in the increase 
of its energy which is found to be more in the case of circular polarized laser. Thus it 
can be inferred that for the same initial laser intensity and static axial magnetic held, the 
circularly polarized laser drives the initially non-resonant particle into cyclotron resonance 
at rate faster than linearly polarized laser. The acceleration of initially non-resonant particle 
caused by its de-phasing from the laser held due to radiative friction is of great signihcance. 
This signihcantly relaxes initial condition in terms of the precise requirement of applied axial 
magnetic held for the resonant acceleration of the particle that has previously hindered the 
experimental realization of the scheme. 

The ehect of radiation reaction force on the characteristic resonance curve is shown in the 
Fig-4. The three energy plots as function of the resonant parameter ‘r' in the hgure corre¬ 
sponds the same normalized static axial magnetic held (.83) and three diherent normalized 
values of the laser amplitudes given by Aq = 10, 300 and 700. From the hgure it can be 
observed that the presence of radiation reaction force broadens the width of the resonance 
curve as well as causes a modihcation in the resonant condition. The resultant resonant 
peaks are shown in the Fig-4 and from the hgure it can be inferred that the width of the 
resonance curve as well as the shift in resonant peaks dehned in terms of ‘r' depends upon 
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the laser amplitude. The shifting in the resonance peaks is similar in its effect to that of 
a damped driven oscillator with a variable coefficient of friction. In the present case the 
oscillator is driven by the axially applied magnetic held together with laser and the dampen¬ 
ing is caused by radiative frictional coefficient which has non-linear dependence upon laser 
amplitude. 

One of the major challenge in the experimental realization of the scheme is the requirement 
of the static magnetic helds of the order of (~ 10^ Tesla) which has been shown to decrease by 
one order (~ 10^ Tesla) by inclusion of radiative friction. In Fig-5, it shown that signihcant 
reduction in magnetic held which are of the order of presently available magnetic helds (~ 10^ 
Tesla) in the laboratory and further gain in particle energy can be simultaneously achieved 
by imparting a small longitudinal momentum (~ Mev) to the initially non-resonant particle. 
From the sub-plots A,B,C and D it can be inferred that for the same initial value of resonant 
parameter ‘r(= .95)', an approximate six times gain in energy 100 TeV) can be achieved 
by four times increase in the initial longitudinal momentum(^ 2.5 — lOMeV). On comparing 
the gain in energy of the particle for linear and circular polarization it is evident that the 
energy gain is more for the circular polarization. 

The initial longitudinal momentum imparted to the particle causes a Doppler down shift 
in the laser frequency as ‘seen' by the particle. As a result the particle remains phase locked 
in the accelerating phase of the laser resulting in large gain in peak energy. In Fig-6, it 
shown that for initially non-resonant particle the peak energy gain depends upon the rate 
of evolution of resonant parameter ‘r'. On comparing hgures (5 and 6), the peak energy 
gain is shown to be more for the slower evolving resonant parameter. The rate of evolution 
of resonant parameter ‘r' is shown to have inverse dependence upon the initial longitudinal 
particle momentum. In Fig-7, it is further shown that the peak energy gain increase linearly 
with the normalized laser amplitude for both linear and circular laser polarization. The 
quantum correction dehned by equation Eq-8 and Eq-9 to the radiation force in the LL 
equation of motion are found to be negligible for various peak energy gains by the particle 
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shown in figure (5 and 6) and corresponding to different parameters. Further a parametric 
space has been identified in this acceleration scheme where the quantum corrected energy 
estimates differ from the one obtained from classical estimates. The effect of quantum 
corrections are shown in Fig-8, which for circularly polarized laser correspond to very high 
laser intensities ( 10^® hF/cm^). 

In this study the energy gain by the particle can be optimized in terms of the laser pa¬ 
rameters i.e its polarization and amplitude along with the initial value of resonant parameter 
‘r', which can be altered by varying particle’s initial longitudinal momentum and applied 
static axial magnetic field. One of the possible experimental ways of extracting the energetic 
particle from laser field is by using a metallic filter which can absorb the laser light and 
whose thickness is much smaller than stopping distance of the particles |39j. 


IV. SUMMARY 

The interaction of charged particle with the combined field of ultra intense laser and 
static axial magnetic field has been explored using Landu-Lifshitz equation. It has been 
shown that the radiation reaction force (or radiative friction) plays a significant role in the 
laser driven Auto-resonant scheme of particle acceleration. The interaction of the stationary 
particle with fields is studied at first, by considering it to be initially in cyclotron resonance 
with applied fields and secondly, by choosing it to be initially non-resonant. For the initially 
resonant case, it has been shown that in contrast to unbounded increase in particle energy 
in the absence of radiation reaction force, the particle energy saturates with the inclusion of 
radiative friction. For the second case corresponding to initially non-resonant particle, the 
inclusion of radiation friction results in the net increase of longitudinal particle momentum 
at each gyration which in turn pushes it to cyclotron resonance resulting in large energy 
gain (~ TeV). The radiative friction eventually causes the de-phasing of the particle from 
the laser field resulting in net saturated increase in the particle energy. It has been shown 
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the presence of radiative friction results in the resonance width broadening as well in the 
shifting of resonance peak to the higher values of the resonant parameter ‘r'. The broadening 
in the resonance width and its shifting is shown to depend upon the laser amplitude. The 
fact that the radiative friction assists in the acceleration of initially non-resonant particle 
is of great significance and is expected to play a role in the experimental realization of the 
scheme. Further, the overall peak energy gain of the particle is shown to have an inverse 
dependence upon the rate of evolution of previously dehned parameter ‘r'. 

In terms of parameter optimization, it has been shown that the peak energy gain of the 
particle for a laser (A ~ Ipm ) with same intensity (~ 10^^can further be enhanced 
(~ 100 TeV) by imparting a small longitudinal momentum (~ 2.5-10 Mev) to the initially 
non-resonant particle prior to the onset of interaction. The initial longitudinal momentum 
causes a Doppler down shift in the laser frequency as ‘seen' by the particle and in turn 
resulting in up to two order reduction in the applied static axial magnetic held requirements 
(~ 200 — 400 Tesla) for the cyclotron resonance and are available in laboratories. It has also 
been demonstrated that the peak energy gain by the particle grows linearly with the laser 
amplitude. The quantum corrections in the energy estimates arising due to the weakening of 
the radiative friction are found to be negligible for the above set of parameters. It has been 
shown that the pronounced quantum corrections to the energy estimates in this acceleration 
scheme will be begin at very high laser intensities (> 10^^IF/cm^). 
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FIGURE CAPTIONS 


Fig(l): The momentum and configuration space particle trajectories for the initially resonant 
in combined field of Circularly Polarized laser (CP) together with static axial magnetic field 
with and without radiation reaction force. 

Fig(2):The momentum and configuration space particle trajectories for the initially non¬ 
resonant in combined field of Circularly Polarized laser (CP) together with static axial mag¬ 
netic field with and without radiation reaction force. 

Fig(3): The energy gain of the initially resonant as well as non-resonant particle in combined 
static axial magnetic field along with that of linearly (LP) and circularly polarized laser (CP) 
respectively. 

Fig(4): Resonant curve for different laser amplititudes. 

Fig(5): Energy gain of the initially non-resonant particle as function of resonant parameter 
‘r' for linearly polarizied (LP) and cricularly polarizied laser (CP). 

Fig(6): Evolution of resonant parameter ‘r' for different initial values of Longitudinal mo¬ 
mentum for lineraly (LP) and circularly polarizied (CP) lasers. 

Fig(7): Peak energy gain as a function of laser amplitude (Aq). 

Fig(8): The comparison of classical and quantum corrected energy estimates. 
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Resonant Momentum and Configuration Particle Trajectory (Polarization:CP) 


Fig:A (nc/o)=1; Ao=50: Pxo=0) Fig:B (nc/(o=1; Ao=50: Pxo=0) 





FIG. 1. The momentum and configuration space particle trajectories for the initially resonant in 
combined field of Circularly Polarized laser (CP) together with static axial magnetic field with and 
without radiation reaction force. 
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Non-Resonant Momentum and Configuration Particle Trajectory (Polarization:CP) 


Fig:A (^(7(0=0.7; Ao=50; Pxo=0) 


Fig:B (nj7o)=0.7; Ao=50; Pxo=0) 



FIG. 2. The momentum and configuration space particle trajectories for the initially non-resonant 
in combined field of Circularly Polarized laser (CP) together with static axial magnetic field with 
and without radiation reaction force. 
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Energy Gain In Auto-Resonant Acceleration With Radiation Reaction Force 
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Fig C (Non-Resonant; £l^m=.7; A(,=50; Pxq=0; LP) 
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FIG. 3. The energy gain of the initially resonant as well as non-resonant particle in combined static 
axial magnetic field along with that of linearly (LP) and circularly polarized laser (CP) respectively. 
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Energy gain as a function of resonant parameter ’r’ 



FIG. 4. Resonant curve for different laser amplititudes. 
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Energy Gain As Function of r(=n^A) 
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Fig :D (Ao=50,Pxo=20,ac/o)=0.023735175) 
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FIG. 5. Energy gain of the initially non-resonant particle as function of resonant parameter V' 
for linearly polarizied (LP) and cricularly polarizied laser (CP). 
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Evolution of resonance parameter r(=n^A) 
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FIG. 6. Evolution of resonant parameter ‘r' for different initial values of Longitudinal momentum 
for lineraly (LP) and circularly polarizied (CP) lasers. 
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Energy Gain As Function of Laser Amplitiude (il(/(a=0.04738184; Pxq=10) 



FIG. 7. Peak energy gain as a function of laser amplitude (^o)- 
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Comaprison between Classical and QED modified RR (Aq=4000) 
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FIG. 8. The comparison of classical and quantum corrected energy estimates. 
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